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Measurements of the angular spectrum of light transmitted through turbid slabs with monodispersions 
of polystyrene spheres have been performed. The results obtained are compared with theoretical cal- 
culations, based on the small-angle approximation of the radiative transfer theory. The experimental 
data and the theoretical results coincide with a high accuracy, which allows us to develop the laser 


diffraction spectroscopy of optically thick light-scattering layers. 
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1. Introduction 


The analysis of the Fraunhofer diffraction patterns of 
different light-scattering media is a standard tool in 
the field of particulate media characterization. 
Both the size? and the shape? of particles can be 
retrieved from angular dependencies of the scattered 
light at small observation angles. For instance, it 
follows for polydispersions of large spherical particles 
(ka >> 1, 2|m — 1|ka > 1, k = 207/), where à is the 
wavelength, m is the relative refractive index, and a 
is the radius of particles),? 


1(6) = Co”? T a*f(a)J(kad)da, (1) 


0 


where C = const, f(a) is the particle size distribution, 
I(8) is the measured scattered intensity, 6 is the scat- 
tering angle, and J,(k0a) is the Bessel function. 
There are many methods to retrieve the particle size 
distribution f(a) from Eq. (1).2 They have already 
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been incorporated into the software of modern parti- 
cle sizers.14-6 

Equation (1) can be applied only in the case of thin 
turbid layers, when the multiple light scattering is 
negligible. However, there is a need for information 
on particle size distributions in dense media with 
strong multiple light scattering. This is because 
some media (e.g., in solid state) cannot be diluted or 
are not under control (e.g., various aerosols and 
clouds in the sky). The inverse problem for multiply 
scattering media was formulated and solved in Refs. 
7-11. 

Our task in this paper is the experimental study of 
the influence of the multiple light scattering on 
Fraunhofer diffraction patterns of polysterene spher- 
ical particles. Experiments were carried out with 
the standard commercially available Sympatec Par- 
ticle Sizer.> Particles were almost monodispersed 
and characterized by mode radii near 5 and 50 pm. 
The results of experiments were compared with the- 
oretical solutions, obtained in Ref. 12. 


2. Experiment 


The experimental setup is presented in Fig. 1. A 
He-Ne laser beam with the wavelength \ = 0.6328 
um passes through a beam expander and a sample 
with monodispersions of polystyrene spheres. The 
diffraction pattern is measured with the help of the 
detector with 31 semicircular elements of a different 
width.” Polystyrene spheres with certified mean di- 
ameters dọ equal to 100 + 2.0 um and 9.685 + 0.064 
um have been used in the experiment. The coeffi- 
cient of variation of the particle size distribution was 
equal to 0.042 for large particles and 0.014 for 
smaller ones. Particles were composed of polysty- 
rene DVB (divinylbenzene, 4—8%) and had the spe- 
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Fig. 1. Experimental setup. 


cific gravity 1.05 g/ml. The refractive index of 
particles at the wavelength 0.59 um was equal to 
1.59. 

The results of the experiment are presented in 
Figs. 2 and 3 for different optical thicknesses t and 
sizes of particles. Diameters of particles were much 
larger than the wavelength of the incident light. 
Note that it follows in a good approximation for the 
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Fig. 2. Angular dependence of the light transmitted through a 

sample with polystyrene spherical particles according to the ex- 

periment and the theory at the optical thickness equal to 0.144, 

0.463, and 3.108. The diameter of spheres is equal to 9.685 um. 
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Fig. 3. Same as Fig. 2 but with diameter of spheres of 100 um and 
optical thickness equal to 0.06, 3.257, and 5.345. 
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extinction efficiency factor Q.., = 2 in this case. 
Thus the optical thickness is approximately given 
by! + = 3Lc/d.,, where where d.p is the Sauter di- 
ameter of particles, L is the geometrical thickness of 
a sample, and c is the volumetric concentration of 
particles in the suspension. It follows for the value 
of c that 


c =a /3L. 


One can see that the value of c is simply expressed by 
the optical thickness of the medium 7. The value of 
t was measured independently during the experi- 
ment. The geometrical thickness L was equal to 20 
mm and d.f ~ do, owing to the high monodispersity of 
particles in a sample. 

The experiment (symbols) shows that the increase 
in the concentration of particles leads to higher val- 
ues of the transmitted intensity in the region of small 
values of t. However, the transmitted intensity re- 
duces with further growth of the optical thickness. 

Indeed, the intensity of the transmitted diffused 
light I4 is equal to zero if there are no particles in a 
cuvette. It grows with the optical thickness ini- 
tially. However, it is equal to zero for semi-infinite 
media. Thus it should be a maximum in the depen- 
dencies J,,(t) at some value oft. This is in correspon- 
dence with our experimental results. It should be 
noted that oscillations in Fig. 3 confirm the high 
monodispersity of particles in a sample. 

Another interesting feature is related to the broad- 
ening of the diffraction pike due to the multiple light- 
scattering phenomenon. The larger the optical 
thickness 1, the larger the half-width of the angular 
spectrum of the transmitted light. It should be 
noted that the half-widths of the small-angle peaks of 
the scattered light also increase with the decrease of 
particle size [see Eq. (1)]. Thus the presence of mul- 
tiple light scattering could lead to underestimation of 
particle size in the standard inversion procedures, 
which do not account for multiple light scattering. 

Figures 2 and 3 differ owing to the absence of min- 
ima in Fig. 2 for low concentrations of particles. 
This is due to the limitation of observation angles by 
4° in the experiment. The first minimum should 
appear at the scattering angle 


_ 3.83 X 180 


d 
TX 


degrees according to the Fraunhofer theory.!4 Here 
x = 2ma/) is the size parameter. The value of 0, is 
approximately equal to 0.4° ata = 50 um. But itis 
~4° ata = 5 um, which is almost outside the mea- 
sured angular range. 

It should be noted that multiple light scattering 
reduces oscillations and makes diffraction patterns 
more smooth. The same effect could be due to the 
polydispersity of particles in a turbid medium. 


3. Theory 


The diffused transmitted intensity can be found from 
the solution of the radiative transfer equation. 1213-15-17 
This equation is written in the following form in the 
case of the normal illumination of a scattering plane— 
parallel layer with randomly oriented particles by a 
plane wave, 


ð T 
os ð Ma =I (7; 0) + | I(t, 3’) p(s, 9’) 
dt 2 


0 


x sin ð'dð’, (2) 


where 9 is the observation angle; wo = Osca/Sext is the 
single-scattering albedo; o,,,, and o,,; are scattering 
and extinction coefficients; T = o,,,l is the optical 
thickness; L is the geometrical thickness of a layer; 
pð, 9’) = % f?" p(®)d@ is the phase function, aver- 
aged on the relative azimuth g; Zis the light intensity; 
and 8 = arccos(cos ð cos 0’ + sin ð sin 8’ cos @) is the 
scattering angle. It should be noted that Eq. (2) does 
not account for the spatial width and the angular 
spread of the incident beam. Thus it can be used 
only if the beam cross section is sufficiently large, 
which is the case for our experiment (see Fig. 1). 

Our primary task here is to study the influence of 
multiple light scattering on Fraunhofer diffraction 
patterns in optically dense disperse layers. Thus we 
consider the behavior of the function I(t, 3) at small 
values ð and assume that 


p Ue 9) _ UG, 8) 
dtr dtr 


cos 


It follows from Eq. (2) in this case that 


dI (7, ð) 
dr 


= -I (r, ð) 


$ zl I(z, 9p, 8')sin 8'd0’. 
0 


(3) 


The integral in Eq. (3) accounts for the multiple scat- 
tering of photons from all scattering directions to the 
direction specified by the observation angle ð. It 
should be noted that the function I(t, 3’) is charac- 
terized by a rather sharp peak at 0’ ~ 0 (see, e.g., the 
results of our experiment in Figs. 2 and 3). The 
same is true for the phase function p(8), which is 
highly peaked in the forward-scattering direction for 
particles large compared with the wavelength of the 
incident radiation. This allows us to change the 
limit of integration in the integral term of Eq. (3) from 
a to © and obtain the following approximate trans- 
port equation, 


dI(t, ð) 
dr 


Sie Hye a | I(r, 8')p(8, 8)8'd9", 
0 


(4) 


where we accounted for the fact that sin 3’ ~ 9’ as 
ð’ — 0 and the phase function p(6) is normalized by 
the following condition: 


1 œ 
al p(6)ede = 1. (5) 


0 


One should expect that approximation (4) provides 
a high accuracy at small values of ð, which are of 
interest in this paper. Equation (4) can be exactly 
solved with several different techniques. For in- 
stance, one can perform the Fourier transform of this 
equation and obtain?2 


In f“ 
I(t, 8) = A exp{—t[1 — wP (o) }}Jo(od)odo, 
Í 6) 


assuming that the intensity of the incident wave is 


where 
1 co 
(ð) = =| J(ad)odo (8) 
QT 3 
is the delta function and I, is the constant. Here 
J (a8) is the Bessel function. The function 
1 o0 
P(o) = l P(9)eJ (08) 0d0 (9) 
0 


in Eq. (6) is the Fourier—Bessel transform of the 
phase function p(6). It is called the space-frequency 
decay function. One can see that P(o = 0) =1lasa 
result of the normalization condition (5). It is easy 
to verify solution (6) by means of substituting it into 
Eq. (4). 

Thus Eq. (6) allows us to investigate diffraction 
patterns of multiply scattering layers by simple 
means. Let us subtract the coherent component of a 
light field [)5(8)exp(—r) [see Eq. (8)] from the general 
solution (6). Then it follows for the diffused trans- 
mitted light J,(7, 5) that 


I ~ 
It, 8) = = exp(—T) | {exp|woP(o)t] — 1} 
0 

x Jo(od)odo. (10) 


This equation can be simplified in the case of thin 
layers. It follows from Eq. (10) as t > 0 that 


I(t, d) = (Io/477) wotp(9), (11) 
where [see Eq. (9)| 


p(s) = 2 T P(o)Ja(oð)odo. (12) 


0 
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Fig. 4. Dependence F(z) at different values of optical thickness 
equal to 0.1, 1, 3, and 5. 


Thus the small-angle diffused intensity is determined 
by the angular distribution of singly scattered light 
as it should be for small values of the optical thick- 
ness. 

Equation (10) allows us to study the dependence of 
diffraction patterns on the optical thickness and con- 
centration of particles. To do so, we should know the 
function P(o) in Eq. (10). That function differs for 
particles of different shapes and internal structures. 
For instance, it follows in the framework of the 
Fraunhofer approximation for monodispersed spher- 
ical particles!® that 


AJ ,°(0x) 


ro (13) 


p(9) = 


Oo = 


2 , 
From Eqs. (9) and (13) one can obtain 


2 o o o \ |?" o 
rotti THE 
(14) 


where u(o/2x) = 1 ato < 2x and u(o/2x) = 0 at ø > 
2x. Thus the diffused intensity [,(7, ®) can be pre- 
sented in the following form [see Eq. (10)]: 


Ie, ð) = DF(r, a), (15) 


where D = 2x°lo/7 is the constant, which does not 
depend on the concentration of particles, and the 
function 


F(t, z) = exp(—7) f (7 [arccos(y) 


-y(1 -ya = 12y (16) 


depends only on the optical thickness 7 and the pa- 
rameter z = 20x. The function F(z) at t = 0.1, 1, 3, 
and 5 is presented in Fig. 4. One can see that the 
intensity of the diffused transmitted light increases 
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Fig. 5. Dependence of transmitted diffused light intensity [see 
Eqs. (15) and (16)] on optical thickness at z = 1.5. 


with the optical thickness, reaches the maximal val- 
ues, and starts to decrease. The amplitude of oscil- 
lations decreases with the optical thickness r. 
However, the half-width of the small-angle peak in- 
creases with the optical thickness. All these fea- 
tures are consistent with the experimental data 
presented in Figs. 2 and 3. 

The dependence of light intensity [Eq. (15)] on the 
optical thickness at z = 1.5 is presented in Fig. 5. 
One can see that the maximum of the intensity at the 
fixed angle occurs at t ~ 1 in this case. It should be 
noted that the function J,(7, 5) in Fig. 5 was normal- 
ized to 1 at t = 0.1. 


4. Comparisons of Computed and Experimental 
Results 


Let us compare results of calculations with the ex- 
perimental data obtained. To this end results of cal- 
culations with Eqs. (15) and (16) were matched to the 
experimental data at 0 = 1° for smaller particles (see 
Fig. 2) and at 6 = 0.1° for larger particles (see Fig. 3). 
The agreement between theory and experiment is 
excellent for all optical thicknesses except tails of the 
angular distributions of transmitted light (see, in 
particular, Fig. 3). Thus simple Eq. (16) can be in- 
deed used for solution of both the inverse and the 
direct problems of light-scattering media optics. 

The poor agreement of the theory and experiment 
at 0 > 2° in Fig. 3 is due to the fact that multiple 
scattering of light refracted and reflected by 
particles!*17 becomes important in this case. This 
is not accounted for in Eq. (16), where the geometrical 
optics component of the scattered light is simply ne- 
glected. Accounting for the geometrical optics scat- 
tering becomes important for larger particles at 
smaller scattering angles than it is for smaller par- 
ticles (see Figs. 2 and 3). This is due to smaller 
half-widths of Fraunhofer diffraction patterns for 
larger scatterers. 

Accounting for the geometrical optics scattering is 
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Fig. 6. Angular dependence of light intensity according to the 
experiment and the theory for samples with polystyrene spheres of 
100-um diameter with accounting for the geometrical optical scat- 
tering at optical thickness equal to 3.257 and 5.345. 


rather simple. Indeed, neglecting the interference 
between diffracted, reflected, and refracted beams, 
one obtains16 


Tsa p70) T Osea 8) 


d g 
Osca R Osca 


p(0) = ; 
where o eat and Caf are parts of the total scattering 
coefficient Csa = Osca? + Osaf, related to diffraction 
and geometrical optics scattering processes, respec- 
tively. 

It follows that oga? = 
particles!8.19 and that 


_ (0) + p“(0) 
a 


The function p?(0) is described by Eq. (13), and the 
function p*(8) represents the geometrical optics con- 
tribution to the phase function p(6). Note that func- 
tions p?(0), p£(0), and p(0) are normalized according 
to Eq. (5). 

Calculations performed in the framework of the 
geometrical optics approximation!>1” show that the 
function p*(6) can be approximated by the following 
simple formula: 


p*(8) = 4B exp(—B6’), (18) 


where the value of B depends only on the refractive 
index of transparent particles and not on their size. 
The multiplier 48 in Eq. (18) is due to the normaliza- 
tion condition (5). 

For the Fourier—Bessel transform (9) of the phase 
function (17) it follows that 


= P%(o) + P8(c) 

2 z 
where the value of P“(c) can be found from Eq. (14) 
and [see Eqs. (9) and (18)], 


P£(o) = exp(—0°/4B). (20) 


Osca for large transparent 


p(9) (17) 
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Results of the comparison of the calculated diffused 
transmitted intensity [see Eqs. (10) and (19)] with 
experimental data for monodipersed spheres with a 
100-um diameter are presented in Fig. 6. One can 
see that accounting for the geometrical optics scat- 
tering [see Eqs. (19) and (20)] indeed provides higher 
accuracy at larger values of 3. Note that the value 
of B = 5 in Eq. (20) was obtained by fitting of exper- 
imental and theoretical results. 

It should be noted that Eq. (10) is also valid for 
randomly oriented nonspherical particles and for 
polydispersed media. However, these media are 
characterized by different spectra P(o) as compared 
with the simple case of monodispersed spheres. 
This results in different angular spectra J,(7, 3) of the 
diffused light. For instance, for disperse media with 
polydispersions of spherical particles it follows that16 


gra) Pléida 
(P(e) = =—- l (21) 
| af (a)da 


0 


where f(a) is the particle size distribution and P(o) is 
described by Eq. (19). 


5. Conclusions 


The experimental measurements of the small-angle 
patterns of particulate media show that the influence 
of multiple light scattering causes the broadening of 
the angular spectrum of the transmitted light and 
reduces the intensity of the Fraunghofer rings for 
monodispersed particles. One can see that Eqs. (15) 
and (16) describe the experimental data with a high 
accuracy. Thus they can be used as a theoretical 
basis for the laser diffraction spectroscopy of turbid 
media. 

The function (P(c)) in Eq. (21) can be retrieved 
from Eq. (10) by the inverse Fourier—Bessel trans- 
form. On the other hand, the Fourier—Bessel trans- 
form of the function (P(c)) [see Eq. (12)] allows us to 
obtain the function p(6) [see Eq. (1)] and reduce the 
inverse problem for a multiply light-scattering me- 
dium to the well-studied case of single scattering. 

It is interesting to note that integral equation (1) 
can be solved analytically by use of the generalized 
Fourier transform.!4 Thus three successive integral 
transforms allow us to obtain the analytical solution 
for the particle size distribution from measurements 
of the diffused transmitted light in the multiple- 
scattering regime. Other possibilities for solving 
the inverse problem associated with Eq. (10) are dis- 
cussed in Ref. 17. 
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